Measurement of Production Properties of Positively Charged Kaons 
in Proton— Carbon Interactions at 31 GeV/c 
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Spectra of positively charged kaons in p+C interactions at 31 GeV/c were measured with the 
NA61/SHINE spectrometer at the CERN SPS. The analysis is based on the fuh set of data collected 
in 2007 with a graphite target with a thickness of 4% of a nuclear interaction length. Interaction 
cross sections and charged pion spectra were already measured using the same set of data. These 
new measurements in combination with the published ones are required to improve predictions of 
the neutrino flux for the T2K long baseline neutrino oscillation experiment in Japan. In particular, 
the knowledge of kaon production is crucial for precisely predicting the intrinsic electron neutrino 
component and the high energy tail of the T2K beam. The results are presented as a function of 
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laboratory momentum in 2 intervals of the laboratory polar angle covering the range from 20 up to 
240 mrad. The kaon spectra are compared with predictions of several hadron production models. 
Using the published pion results and the new kaon data, the K'^ /■k'^ ratios are computed. 

PACS numbers: 13.85.Lg,13.85.Hd,13.85.Ni 
Keywords: p+C interaction, inclusive kaon spectra 



I. INTRODUCTION 

The SPS Heavy Ion and Neutrino Experiment 
(NA61/SHINE) at CERN pursues a rich physics pro- 
gram [1~4]. Hadron production measurements in p+C 
and TT+C interactions will improve calculations of neu- 
trino fluxes in the T2K experiment [5], and simulations 
of cosmic-ray air showers in the Pierre Auger and KAS- 
CADE experiments [6, 7]. The heavy ion program in- 
vestigates p-l-p, p+Pb and nucleus -I- nucleus collisions 
at SPS energies, to study the onset of deconfinement 
and search for the critical point of strongly interacting 
matter. Charged pion spectra in p+C interactions at 
31 GeV/c were already published [8] and used for neu- 
trino flux prediction in T2K [9]. This article presents 
new measurements of positively charged kaon spectra in 
p+C interactions at 31 GeV/c, based on the data col- 
lected during the first running period in 2007. A detailed 
description of the experimental apparatus and analysis 
techniques can be found in [8] . 

T2K - the long baseline neutrino experiment from 
J-PARC in Tokai to Kamioka (Japan) - aims to pre- 
cisely measure the — v^. appearance and disap- 
pearance [5, 9]. The neutrino beam is generated by the 
J-PARC high intensity 30 GeV (kinetic energy) proton 
beam interacting in a 90 cm long graphite target to pro- 
duce TT and K mesons, which decay emitting neutrinos. 
The resulting neutrino beam is aimed towards a near de- 
tector complex, 280 m from the target, and to the Super- 
Kamiokande (SK) far detector located 295 km away at 
2.5 degrees off-axis from the hadron beam. Neutrino os- 
cillations are probed by comparing the neutrino event 
rates measured in SK to the predictions of a Monte-Carlo 
simulation based on flux calculations and near detector 
measurements. Until the NA61/SHINE data were avail- 
able, these flux calculations were based on hadron pro- 
duction models tuned to sparse available data, resulting 
in systematic uncertainties which are large and difficult 
to evaluate. Direct measurement of particle production 
rates in p+C interactions allows more precise and reli- 
able estimates. Presently, the T2K neutrino beam-line 
is set up to focus positively charged hadrons, to produce 
a beam. While charged pions generate most of the 
low energy neutrinos, positively charged kaons generate 
the high energy tail of the T2K beam, and contribute 
substantially to the intrinsic Vf. component in the T2K 
beam. 

Positively charged kaons whose daughter neutrinos 
pass through the SK detector constitute the kinematic 
region of interest, shown in Fig. 1 in the kinematic vari- 
ables p and 9 - the momentum and polar angle of par- 



ticles in the laboratory frame. The low statistics avail- 
able in the 2007 pilot data set imposes a p, 6* binning 
which covers only the most populated region of phase 
space relevant for T2K. Moreover, the statistics of the 
2007 data does not allow for measurements of negatively 
charged kaons. An order of magnitude larger data set 
was recorded in 2009, and, when analyzed, will lead to 
essentially full coverage. The NA61/SHINE data on kaon 
production will allow also to test and improve existing 
hadron production models in an energy region which is 
not well constrained by measurements at present. Sev- 
eral production measurements in this energy range 
were performed previously [10-16]. 




p [GeV/c] 

Figure 1: (Color online) The prediction from the T2K beam 
simulation; the {p, 9} distribution for positively charged kaons 
weighted by the probability that their decay produces a neu- 
trino passing through the SK detector. The binning used 
in the present analysis is superimposed; the kinematic range 
considered is 1.6<p<7.2 GeV/c and 20< 61 <240 mrad. 



II. THE EXPERIMENTAL SETUP 

The NA61/SHINE apparatus (for details see Ref. [8]- 
Sec. H) is a wide acceptance spectrometer made of four 
large volume Time Projection Chambers (TPCs): two 
Vertex TPCs (VTPC-1 and VTPC-2) placed in the mag- 
netic field produced by two superconducting dipole mag- 
nets and two Main-TPCs (MTPC) located downstream 
symmetrically with respect to the beam line (Fig. 2). In 
the forward region, the setup is complemented by a time- 
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of-flight (ToF-F) detector array horizontally segmented 
into 64 scintillator bars read out at both ends by photo- 
multipliers. The time resolution of each scintillator is 
about 115 ps [17]. For the study presented here the mag- 
netic field was set to a bending power of 1.14 Tm. This 
leads to a momentum resolution a{p)lp'^ in the track 
reconstruction of about 5 x 10"'^ (GeV/c)~^. The tar- 
get used is an isotropic graphite sample with a thickness 
along the beam axis of 2 cm, equivalent to about 4% of 
a nuclear interaction length, Aj. During the data taking 
the target was placed 80 cm upstream of the VTPC-1. 

III. ANALYSIS TECHNIQUE 

This section presents details on data selection and bin- 
ning, the kaon identification method as well as correc- 
tions and systematic errors. 

A. Data binning 

The adopted binning scheme is mainly driven by the 
available statistics and is presented in Fig. 1. Details can 
be found in Table I. The highest 9 limit is determined 
by the requirement for the track to be in the geometrical 
acceptance of the ToF-F detector. 

Only two angular intervals are defined. The lowest 9 
value is set to 20 mrad in order to exclude tracks passing 
close to the edges of the TPCs where the reconstruction 
efficiency is lower and the calculation of the correction for 
the acceptance is less reliable. The first angular interval 
extends up to 140 mrad so as to cover most of the T2K 
relevant 9 range and, combined with a 0.8 GeV/c momen- 
tum bin width, to have a few thousands of selected tracks 
per interval (Table I). Measurements were performed up 
to 7.2 GeV/c. This choice comes from the fact that in the 
relativistic rise region (above 4-5 GeV/c) particle iden- 
tification requires extracting the rapidly decreasing kaon 
signal from the predominant proton one. With the avail- 
able statistics of 2007 data the applied procedure turned 
out to be robust only up to about 7 GeV/c. 

B. Event and track selection 

This analysis is based on 452 x 10'' reconstructed events 
collected during the 2007 data taking. Only events for 
which a beam track is properly reconstructed are se- 
lected. The beam trajectory is measured with a set of 
Beam Position Detectors (BPD) placed upstream of the 
target ([8]-Sec. III,V). Several criteria were applied to se- 
lect well-measured positively charged tracks in the TPCs 
and ensure high reconstruction efficiency as well as to 
reduce the contamination of tracks from secondary inter- 
actions: 

(i) track momentum fit at the interaction vertex should 
have converged. 



(ii) a minimum of 12 reconstructed points in the two 
TPCs used for momentum measurement, VTPC-1 
and VTPC-2, is required, 

(iii) the distance of closest approach of the fitted track 
to the interaction point (impact parameter) is re- 
quired to be smaller than 4 cm in both transverse 
directions, 

(iv) the track must leave the primary vertex at an az- 
imuthal angle (j) within ±20° around the horizontal 
plane, for the first 9 interval, and ±10° for the sec- 
ond; this excludes most of the tracks traversing the 
detector in the regions where the reconstruction ca- 
pability is limited by the magnet aperture or by the 
presence of uninstrumented regions in the VTPCs, 

(v) the track must have an associated ToF-F hit. 

The position of a ToF-F hit is determined only in the x 
direction and with a precision given by the width of the 
scintillator slat producing the signal (~ 10 cm) (Fig. 2). 
A ToF-F hit is then associated to a track if the trajectory 
can be extrapolated to the pertaining slat. 




VTPC-1 VTPC-2 




MTPC 

13.5 m 



Figure 2: (Color online) Schematic drawing of the experimen- 
tal apparatus (Sec. II) and a reconstructed p-l-C interaction 
event. Yellow (green) points indicate TPC clusters (not) as- 
sociated to reconstructed primary tracks. Stars correspond 
to hits reconstructed in the ToF-F. Red lines are the fitted 
particle trajectories. 

The relatively short proper decay length of kaons (cr ~ 
3.7 m) compared to the longitudinal extension of the de- 
tector (~ 13 m) implies that almost 60% of these particles 
produced in the target at the lowest laboratory momen- 
tum considered, 1.6 GeV/c, decay before reaching the 
ToF-F array. 

Thanks to the high Q values of the kaon decay chan- 
nels, kink topologies are correctly reconstructed as a pri- 
mary and a secondary track (i.e. not fitted to the primary 
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vertex) with an efficiency fiiglier tlian 98%. Nevertlieless, 
a ToF-F hit would stiU be associated to a kaon decaying 
in flight if the secondary particle produces a hit along the 
same scintillator bar in which a hit from the primary is 
expected. As a consequence the time of flight measure- 
ment could be significantly biased. Such special topolo- 
gies can be effectively rejected by considering only tracks 
reconstructed until the downstream edge of the MTPCs. 
(see Sec. HID). This can be achieved by the following 
cut: 

(vi) the z position of the last reconstructed point, z/astj 
must differ by less then 50 cm from the position of 
the last potential point (see also Fig. 2). 

The sample can be divided into tracks emitted in the 
same direction in which they are bent by the magnetic 
field (i.e. > for positive charges) and tracks emitted 
in the opposite direction (see Fig. 2). These two inde- 
pendent sub-samples were treated separately. In fact, for 
the same {p, 0} bin the fraction of accepted tracks can 
differ significantly between the two topologies. In order 
to improve the accuracy of the Monte Carlo correction 
and to reduce the dependence on the model used for 
simulation (Sec. HID), only one track topology per bin, 
the one with the highest acceptance, is chosen. 



C. Combined particle identification 

The specific ionization energy loss, dE/dx, measured 
in the TPCs, and the time-of- flight, measured with the 
ToF-F, are used to identify charged particles of different 
masses in the chosen momentum intervals. 

For each track, the dE/dx is defined as the mean cal- 
culated from the lowest 60% of the cluster charges. This 
reduces the effect of the long Landau tail in the single 
cluster charge distribution to the extent that the dE/dx 
experimental resolution can be assumed to be Gaussian 
within the precision which can be achieved with the avail- 
able statistics. 

A dedicated calibration procedure is applied to the 
data to correct the measured cluster charge deposits for 
various detector effects, e.g. charge absorption along the 
drift path, effective sample length and variations of gain 
between the TPC sectors. After the calibration the cor- 
rected data are fitted to a function parameterizing the 
Bethe-Bloch relation [18]. This function can then be used 
to calculate the mean dE/dx (the mean specific energy 
loss) for a given momentum and for a given mass hypoth- 
esis. 

Kaon identification based only on the dE/dx informa- 
tion is impossible in the momentum bins below about 
4 GeV/c, since protons and kaons have similar values of 
dE/dx in this kinematic range (see Fig. 3). 
Another approach is to identify particles through the 
measurement of their mass squared m?. This method 
is based on the measurement of the time-of-flight which 



along with the track length and momentum allows the 
calculation of the particle mass. 

Identification based only on the information is 
effective at low momenta, below about 4 GeV/c, where 
the separation power is higher (Fig. 3) than for dE/dx. 
The uncertainty on the mass measurement is in fact 
dominated by the time resolution, meaning that the 
m? resolution worsens quadratically with increasing 
momentum (Fig. 5). 

Given the complementarity of the two approaches the 
particle identification (FID) capability can be improved, 
over the whole kinematic range, by using the combined 
dE/dx and m? information. 

Figures 3 and 4 show distributions of measured dE/dx 
versus m? in the first angular interval, 20 < 6* < 140 mrad 
for several momentum bins. Accumulations correspond- 
ing to the different particle types are clearly observable 
(in particular for kaons) in contrast to the individual 
dE/dx and distributions. 

The dE/dx — rm? function used to fit the yields of pro- 
tons (p), kaons (i^"*"), pions (Tr"*") and positrons (e+) to 
the data is a superposition of four bi-dimensional Gaus- 
sians. 

A dedicated Monte Carlo simulation was used to vali- 
date the accuracy of the simple Gaussian approximation 
so as to ensure that no biases related to the definition of 
the particle distribution function are introduced. The se- 
lected function has 20 parameters: 4 yield-, 8 width- and 
8 mean-parameters. In all analyzed bins the accumu- 
lations in the dE/dx-rn^ plane are fully separated from 
other particles (Fig. 3) . Thus the number of parameters 
relevant for the kaon yield determination reduces to 15. 
Arguments and details concerning the selected function 
and further reduction of the number of fitted parameters 
are briefly outlined below. 

(i) Width oj the dE/dx distributions. 

The measured width of the different dE/dx peaks is 
the result of the experimental (detector) resolution 
plus the smearing due to the variation of the dE/dx 
mean value with the momentum. The second con- 
tribution depends on particle type (mass) and is 
non-negligible in the lower momentum bins, where 
the relativistic rise of the energy loss is steep. In or- 
der to separate the two contributions the variance 
of each Gaussian {i = p, K, tt) was expressed as 
the sum of two terms: a? = a'^^p + af . The first 
is the experimental resolution and is treated as a 
free fit parameter. The second, a constant term, 
was derived as the variance of the distribution of 
the dE/dx values calculated from the Bethe-Bloch 
function for all the tracks in the respective {p. 0} 
bin. Despite the fact that momentum distributions 
of protons, kaons and pions differ somewhat, the 
momentum distribution for all particles was used 
to evaluate cr^fcj„. This procedure was justified by 
Monte Carlo studies which showed that, in the con- 
sidered kinematic range, a particle-type-dependent 
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Figure 3: Scatter plots of dE/dx versus measured for the selected tracks in three {p, 6} bins. 
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Figure 4: Example of a bi-dimensional fit to the dE/dx-m? distribution. The m? and dE/dx projections are also shown 
superimposed with the results of the fitted functions. 
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The assumption of a resolution Uexp independent of 
the particle type (mass) was validated by investi- 
gating the dependence of the dE/dx resolution on 
the mean dE/dx. From a fit performed to the pion 
distribution in bins where the accumulation is well 
isolated, no indication of a dependence of the peak 
width on the peak position was found. This result 
allows to use the same experimental resolution for 
all hadron species. 

The experimental resolution is inversely propor- 
tional to the square root of the number of dE/dx 
measurements for a track. The latter vary from 
track to track. Therefore the actual distribution of 
the energy loss is rather a superposition of Gaus- 
sians of different widths. The simple Gaussian ap- 
proximation is still applicable since the topology of 
the selected tracks results in a narrow distribution 
of the number of points peaked above 100 in each 
{p, 9} interval. This finding was validated through 



a dedicated Monte Carlo simulation. 

(ii) Mean of the dE/dx distributions 

Since the mean energy loss depends only on the 
momentum to mass ratio, the dE/dx distributions 
of different particle species with the same momen- 
tum distribution are shifted. In the fit only the 
pion mean energy loss was a free parameter. The 
kaon and proton mean-paramters were instead cal- 
culated using the fitted pion mean and the shifts 
calculated from the Bethc-Bloch parametrization. 

(iii) Width of the m? distributions 

The Gaussian approximation is adopted also for the 
projected m? distributions. 

The MC simulation shows that the distortion of the 
Gaussian shape due to the increase of the vari- 
ance with the momentum has a negligible impact 
on the fitted yields. The m? width-parameters of 
p's, K's and tt's are fitted independently. 

(iv) Mean of the m? distributions 

The mean-parameters of p's, JC's and tt's are 
fitted independently. 
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Figure 5: The resolution versus the measured momentum 
for tracks with 20< 6 <140 mrad. Parabolic fits to the best 
measured points are superimposed. 
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Figure 6: Bethe-Bloch function fit to mean dE/dx values 
(data points). The points used in the fit are from protons 
and kaons of < 6 and pions of > 20. 



Thus the fitted parameters relevant for the kaon yield 
determination are: three yield-parameters as well as one 
width- and one mean-parameter for the dE/dx Gaus- 
sians, as well as three width- and three mean-parameters 
for the Gaussians. 

Performance studies carried out on simulated data 
show that, in the higher momentum bins where the 
K/p and K/tt production ratios rapidly decrease and the 
peaks significantly overlap, the results depend strongly 
on the details of the fit procedure, e.g. initial values and 
bounds of the fitted parameters. This in particular con- 
cerns the dE/dx and mean-parameters as well as the 
vn? width-parameters. Along with the definition of the 
fit function, it is therefore mandatory to establish a pro- 
cedure for the precise determination of the initial values 
and bounds of the fitted parameters. 



For the distributions in m?, the strategy is to esti- 
mate mean- and width-parameters at high momenta by 
extrapolating the fit results from low momentum bins. 

First, mean-parameters of protons, kaons and pions are 
fitted in the bins in which particle accumulations are well 
separated (e.g. up to 4 GeV/c for 20<6'<140 mrad). The 
fitted values are independent of the momentum, which 
proves the correctness of the ToF-F calibration proce- 
dure. The average of the fitted mean-parameters and its 
error are used as the initial value and bounds, respec- 
tively, in the final fit. Second, the m? width-parameters 
are fitted over the whole momentum range. Figure 5 
shows the results for protons and pions. As expected, 
the momentum dependence of the m? width-parameter 
is well fitted by a second order polynomial. The width- 
parameter is larger for protons than pions, since for the 
same momentum and path length, the time of flight is 
larger for heavier particles. The kaon width-parameter 
cannot be fitted because of the small separation in m? of 
the kaon and pion Gaussians. It is therefore calculated 
using the fitted pion width at the momentum bin cen- 
ter rescaled for the correct mass value. The proton and 
pion resolution functions are derived via parabolic fits ex- 
tended over the momentum range where peaks are well 
separated (e.g. up to 5.6 GeV/c for 20< 9 <140 mrad, 
see Fig. 5). Initial values for width-parameters are there- 
fore calculated by extrapolating the resolution functions 
to the center of the momentum bin. Extrapolation errors 
are calculated as well and used to constrain the allowed 
fit parameter range. 

Figures 3 and 4 show that in the dE/dx-w? plane the 
peak of the pion accumulation can be well determined in 
all analyzed bins. For protons and kaons this is possible 
in several bins. Thus a precise fit of the dE/dx mean- 
parameter is possible over the whole range relevant for 
the analysis. The results are presented in Fig. 6. Points 
with low and intermediate /?7 correspond to protons and 
kaons while points with high correspond to pions. 
Thus, points measured with high precision are the ones 
with lower and higher j3"f while the intermediate region 
corresponds to protons and kaons of high momenta for 
which the accumulations overlap (Fig. 3). 

This means that, for intermediate (3^ values, an ac- 
curate parametrization of the expected dE/dx can be 
obtained from a Bethe-Bloch curve fitted only to the 
high precision measured points, which are in the low- 
and high- (3^ region, as specified in Fig. 6. 

The precision of this procedure, estimated from the 
residuals with respect to the fit points, turned out to 
be about 0.5%. This value meets the requirement, de- 
rived from the Monte Carlo studies, of a precision be- 
low 1% needed to keep the systematic error below the 
statistical uncertainty. Estimates of the systematic er- 
rors are reported in Sec. HIE. It is worth noting that 
the required precision on the determination of the dE/dx 
mean-parameters cannot be achieved by using the Bethe- 
Bloch curve as calculated from the calibration procedure 
since this represents an overall parametrization of data 
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averaged over the detector as a whole covering a wider 
range of track topologies. 

Finally, an example of the function fitted to data using 
the Maximum Likelihood method is shown in Fig. 4. 

D. Correction factors 

The Monte Carlo simulation described in Refs. [8]-Sec. 
IVC and [19] was used to calculate corrections for kaon 
decay, secondary interactions in the target and detector 
material and track reconstruction efficiency. Two strate- 
gies were implemented to correct the raw spectra: 

(i) the different biases were calculated separately with 
the MC and applied successively to the data, 

(ii) a global Monte Carlo correction taking into account 
all the above effects. 

The MC correction factors calculated with the first strat- 
egy are shown in Fig. 7 for the first angular bin. The 
separation of corrections clearly shows that the decay 
correction dominates. This correction was estimated, for 
each {p, 9} bin, by computing the fraction of kaons pro- 
duced in the primary interaction which reach the ToF-F 
wall before decaying. As such, this procedure assumes 
that the track selection defined in Sec. IIIB, in particu- 
lar the cut on the z position of the last measured point 
(ziast), is fully efficient in the isolation of a pure sam- 
ple of stable kaons. Before applying this correction it 
is therefore necessary to correct the sample of identified 
kaons for the contamination from those kaons which de- 
cayed before reaching the ToF-F but were still associated 
to a ToF-F hit (see Sec. IIIB item (v) and (vi)). This 
contamination was found to be of about 2-3%. 

The feed down correction concerns kaon tracks fitted 
to the primary vertex but not produced in the primary 
interaction. Feed down is only observed from kaons pro- 
duced by secondary interactions in the target and the 
effect is below 2% over the whole kinematic range. 

The ToF-F efficiency was accounted for and calculated 
from the data by requiring that a track traversing the 
ToF-F wall generates a hit in the ToF-F. 

Except for the first momentum bin, the small correc- 
tion for geometrical acceptance (on top of the factor due 
to the selection of the 4> wedge) reflects the fact that 
only full acceptance regions were selected as described 
in Sec. IIIB. Tracks are also required to hit ToF-F. For 
these long tracks the reconstruction efficiency is close to 
unity. 

The contribution related to kaon losses due to sec- 
ondary interactions in the detector material is about 2%. 

Within the second strategy, global correction factors 
were calculated by means of different detector response 
unfolding methods [20]. All the applied methods re- 
turned results whose discrepancies are significantly below 
the associated statistical errors. 
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Figure 7: (Color online) Monte Carlo correction factors as 
a function of momentum calculated for the angular interval 
20 < 6* < 140 mrad. 

Correction factors obtained through the unfolding 
technique were compared to those calculated as the prod- 
uct of the individual factors derived from method (i). 
Results turned out to be compatible within the statisti- 
cal fluctuations and no systematic biases were observed. 
This implies that all the relevant contributions to the cor- 
rections were identified and quantified correctly and that 
the different contributions are weakly correlated. More- 
over, since the unfolding procedure accounts for event 
migration between adjacent bins, the compatibility of the 
two approaches proves that bin-by-bin correlation of the 
correction factors is negligible as well. 

E. Systematic errors 

The main contribution to the overall systematic error 
on the yield correction arises from the decay correction. 

Three main contributions to the systematic uncer- 
tainty of the decay correction are discussed below. 

(i) The contribution due to kaons decaying in flight. 
Most of the unstable kaons in the selected final 
sample are decays which take place in the region 
of space delimited by ziast and the ToF-F surface. 
An alternative way to derive the decay correction is 
then to calculate the survival probability only until 
z = ziast- In this procedure one must correct for 
the small contamination of kaons decaying between 
Ziast and the ToF-F which are not associated to a 
ToF-F hit. The two methods give results which dif- 
fer by less than 1%. The systematic error assigned 
to the bias related to kaon decays was therefore con- 
servatively calculated as 50% of the corresponding 
contamination (see Sec. HID). 

(ii) The contribution due to the model dependence of 
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momentum distributions. 

The value of the decay correction depends on the 
momentum distribution within each {p, 9} bin and 
therefore on the specific MC model used for the 
event generation. Significant discrepancies may be 
observed especially at low momenta where the de- 
cay probability (the correction) is larger and the 
momentum distribution is steeper. This effect was 
quantified by using various MC models and com- 
paring the resulting correction factors with those 
calculated directly from the data. More precisely, 
the latter factor was derived for the three low- 
est momentum bins, where unique identification of 
kaons is possible on a track-by-track basis: the raw 
yield was corrected by re-weighting each track with 
its decay probability, calculated for the measured 
momentum and track length. The resulting devia- 
tion of the correction factor was negligible except 
in the first momentum bin where a maximum dif- 
ference of about 2% was found. 

(iii) The contribution due to the uncertainty in the re- 
constructed track length. 

Since only well measured tracks which traverse the 
entire spectrometer are retained, a precision of a 
few millimeters is achieved on both track length 
calculation and extrapolation to the ToF-F surface. 
This translates into a negligible (order of 0.1 %) er- 
ror on the decay correction factor. 

Other sources of systematic errors include uncertain- 
ties on the ToF-F and reconstruction efficiency and on 
the contribution from secondary interactions. 

Systematic uncertainties on the ToF-F efficiency come 
from the accuracy of the calibration procedure. The es- 
timated value is 2%. 

The systematic error on the reconstruction efficiency 
was estimated by varying the track selection cuts. The 
induced bias is small compared to the statistical fluctua- 
tions; this is expected considering the high reconstruction 
efflciency. 

Finally a systematic error corresponding to 30% of the 
correction was assigned to the contribution of secondary 
interactions [8]. 

Systematic uncertainties related to the procedure used 
for the FID (Sec. Ill C) were quantified by studying the 
dependence of the fitted kaon yields on the input fit pa- 
rameters, namely: the central values of the dE/dx and 
rn^ distributions and the widths in m^. The systematic 
error was obtained by varying both the input parameters 
and limits and considering the subsequent variation of 
the returned particle yield. 

In particular, the relative distances of the proton, kaon 
and pion dE/dx peaks were varied by an amount cor- 
responding to the largest of the residuals between the 
fitted Bethe-Bloch curve and the measured points (see 
Sec. IIIC). The allowed range for the peak posi- 
tion and resolution were enlarged by a factor of three 
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Figure 8: (Color online) Breakdown of systematic er- 
rors as a function of momentum for the angular interval 
20 < 6* < 140 mrad. 

with respect to the original value which had been de- 
rived from the errors in the parametrization as described 
in Sec. IIIC. Results are shown in Fig. 8. 

As underlined in the introduction of the Sec. Ill, the 
sensitivity to the kaon signal decreases with the momen- 
tum, therefore, the fitted kaon yield depends significantly 
on the definition of the input parameters. This explains 
the steep increase of the systematic error in the last mo- 
mentum bin. 

The relative contributions of all considered systematic 
errors are shown in Fig. 8 for the first angular interval. 



IV. RESULTS 

The K'^ spectra presented in this paper refer to posi- 
tively charged kaons produced in strong and electromag- 
netic processes in p+C interactions at 31 GeV/c. Differ- 
ential inclusive cross sections were derived following 
the procedure described in [8, 21]. 

The results are presented in Fig. 9 and Table I as a 
function of momentum in the two considered intervals of 
polar angle. Momentum and polar angle are calculated 
in the laboratory system. 

The ratio of to 7r+ production cross sections is 
shown in Fig. 10. The 7r+ spectra are taken from Ref. [8], 
central values and errors are recalculated to match the 
binning used in the current analysis. Numerical values 
are reported in Table I. 

V. COMPARISON WITH MODEL 
PREDICTIONS 

In this section we compare the spectra in p-l-C in- 
teractions at 31 GeV/c with the predictions of hadronic 
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Table I: The NA61/SHINE results for the differential production cross section in the laboratory system, da /dp, and the 
K+ to TT^ ratio of production cross sections, /cr'^ , for p+C interactions at 31 GeV/c. Each row refers to a different 
[piow < P < Pup, dlow < ^ < Sup) bin, where p and 9 are the kaon momentum and polar angle in the laboratory frame. is the 
total number of selected tracks and Nk is the fitted raw number of kaons. The central value as well as the statistical (Astat) 
and systematic (Asys) errors of the cross section are given. The overall uncertainty (2.5%) due to the normalization procedure 
is not included. For the if^/vr^ ratio, errors are calculated taking into account only statistical uncertainties. 



> 
a 



3 



1- 



20<e<140 mrad 



6 

p [GeV/c] 



140<e<240 mrad 



p [GeV/c] 




6 

p [GeV/c] 



Figure 9: Differential cross sections for K'^ production in 
p+C interactions at 31 GeV/c. The spectra are presented 
as a function of laboratory momentum, p, in two intervals of 
polar angle, 6. Error bars indicate statistical and systematic 
uncertainties added in quadrature. The overall uncertainty 
(2.5%) due to the normalization procedure is not included. 



Figure 10: Ratio of K'^ to 7r+ production cross sections in 
p+C interactions at 31 GeV/c. The ratios are presented as a 
function of laboratory momentum, p, in two intervals of polar 
angle 9. Errors are calculated taking into account only statis- 
tical uncertainties. Predictions of hadron production models 
are superimposed. 



event generators. Three models, Venus4.12 [22, 23], 
FLUKA2008 [24], and UrQMDl.3.1 [25, 26] were selected 
for this purpose. They are part of the CORSIKA [27] 
framework and are commonly used for the simulation of 
hadronic interactions at energies below 80 GeV in exten- 
sive air showers [28]. Venus is also the standard model 
for acceptance simulations of the NA49 and NA61 Col- 
laborations. 

These models were already compared to the pion spec- 
tra measured by the NA61/SHINE Collaboration [8]. 
Motivated by this comparison, a correction of a technical 



shortcoming of the UrQMD model was proposed in [26] . 
This correction does not affect the kaon spectra presented 
here. Therefore we compare the data with the original 
implementation. 

The results of the comparison between data and mod- 
els are presented in Fig. 11. In order to avoid un- 
certainties related to the different treatment of quasi- 
elastic interaction and to the absence of predictions for 
inclusive cross sections, spectra are normalized to the 
mean multiplicity in all production interactions. For 
the data, the normalization relies on the p+C inclusive 
production cross section aprod which was found to be 
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229.3±1.9±9.0 mb. The production cross section is cal- 
culated from the inelastic cross section by subtracting 
the quasi-elastic contribution. Therefore production pro- 
cesses are defined as those in which only new hadrons 
are present in the final state. Details of the cross section 
analysis procedure can be found in [8, 21]. 

The qualitative behaviour of the data is well repro- 
duced by all models. The quantitative differences can be 
related to the two main production processes of kaons: 
pairwise production of a together with another K 
meson and production of a K'^ together with a A baryon. 
The latter process dominates kaon production at large 
momenta and small angles due to the leading particle ef- 
fect. K'^ at large angles and low momenta stem from pair 
production of K mesons. Both Fluka and Venus pro- 
vide a reasonable description of the pair-produced kaons. 
On the other hand, none of the models is in full agreement 
with the small-angle data. While the Venus model over- 
estimates the production of at small angles, Fluka 
and UrQMD predict a slightly lower kaon production rate. 

Figure 10 also shows the comparison between models 
and data for the ratio of to 7r+ production cross 
sections: UrQMD is in good agreement with the data, 
Fluka provides a reasonable description, while Venus 
overestimates the production cross section ratio for both 
small and large angle intervals. 
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Figure 11: Comparison of measured spectra with model 
predictions. Distributions are normalized to the mean 
multiplicity in all production p-f C interactions. The vertical 
error bars on the data points show the total (stat. and syst.) 
uncertainty. The horizontal error bars indicate the bin size in 
momentum. 



VI. SUMMARY 

This work presents measurements of differential pro- 
duction cross sections of positively charged kaons in p-l-C 



interactions at 31 GeV/c. These data are essential for 
precise predictions of the high energy tail and intrin- 
sic electron neutrino component of the initial neutrino 
fiux for the T2K long baseline neutrino oscillation exper- 
iment in Japan. Furthermore, they provide important in- 
put to improve hadron production models needed for the 
interpretation of air showers initiated by ultra high en- 
ergy cosmic particles. The measurements were performed 
with the large acceptance NA61/SHINE spectrometer at 
the CERN SPS. A set of data collected with a 4% Ai 
isotropic graphite target during the pilot NA61/SHINE 
run in 2007 was used for the analysis. Positively charged 
kaon spectra as a function of laboratory momentum in 
two intervals of the polar angle were obtained. The final 
spectra were compared with predictions of hadron pro- 
duction models. The data presented in this paper are 
already provided to T2K for the calculation of the neu- 
trino flux. Meanwhile, a much larger data set with both 
the thin (4% Ai) and the T2K replica carbon targets was 
recorded in 2009 and 2010 and is presently being ana- 
lyzed. 
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